On the feasibility of studying vortex noise in 2D superconductors with cold atoms 
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We investigate the feasibility of using ultracold neutral atoms trapped near a thin superconductor 
to study vortex noise close to the Kosterlitz-Thouless-Berezinskii transition temperature. Alkali 
atoms such as rubidium probe the magnetic field produced by the vortices. We show that the 
relaxation time Ti of the Zeeman sublevel populations can be conveniently adjusted to provide long 
observation times. We also show that the transverse relaxation times T2 for Zeeman coherences are 
ideal for studying the vortex noise. We briefly consider the motion of atom clouds held close to the 
surface as a method for monitoring the vortex motion. 
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Ultracold neutral atoms trapped and manipulated on 
atom chips [l| can be used as sensitive probes for a wide 
range of phenomena. Examples include accurate mea- 
surements of gravity Q, imaging the magnetic or elec- 
tric landscape near wires and magnetic films 0, 0, 
measuring the field due to Johnson noise near metal- 
lic and dielectric surfaces [1, 0, H, [1] and probing the 
Casimir-Polder force as a function of atom-surface dis- 
tance [10(. Until recently, the atom chips used in these 
kinds of measurements have been at room temperature 
and the corresponding thermal fluctuations of the mag- 
netic field have caused mixing of the Zeeman sublevels. 
Depending on the material of the surface and its distance 
from the atoms, this spin relaxation time is typically of 
order 1-100 seconds. 

The use of superconducting films at cryogenic tempera- 
tures has been proposed as a way to reduce thermal noise 
flll | and cold atoms have now been trapped near a super- 
conducting surface [l2j]. In theoretical studies of atom- 
superconductor interactions, the arguments employed so 
far have been based on a simple two-fluid model [l3[ 
or on data extracted from surface impedance measure- 
ments on bulk superconductors [llj. These approaches 
apply to three-dimensional superconducting materials, 
but the thin films normally used on atom chips are typi- 
cally closer to two-dimensional objects for which 3D the- 
ory is not strictly valid. Near thin film superconductors, 
the magnetic noise is generated primarily by vortex mo- 
tion, which is absent in 3D superconductors. This raises 
the possibility that atoms trapped near the surface of a 
superconducting atom chip might be able to probe the 
physics of vortices. In this article, we start to explore 
the feasibility of using cold atoms to investigate vortex 
noise in a 2D superconductor. 

Rather than going into the theory of 2D supercon- 
ductors and their vortex dynamics, it is enough for our 
present purpose to be guided by existing experimental 
data. In the experiment reported in |l4|, a SQUID 



loop roughly 1mm in diameter measured the flux noise 
spectrum produced by vortices in a 2D superconduct- 
ing Josephson junction array some 100/zm away. The 
detected flux provided statistical information on the dy- 
namics of the vortices in the Josephson array. Here we 
explore how the observed vortex behaviour would influ- 
ence the magnetic sublevel populations and spin coher- 
ence of ultracold atoms trapped above such a surface. 

Our first goal is to relate the flux noise power mea- 
sured in [l4[ to the expected spin flip rate for a magnetic 
atom trapped near the surface. To this end, we write the 
spectral density of the flux noise (per Hz of bandwidth) 
at angular frequency w as 



S (.f) = 27r j d 2 xd 2 y(B z (x,z;cu)Bi(y, 



z;u)) (1) 



where A is the area of the pick-up loop placed parallel 
to the surface at height z. The points (x, z) and (y, z) 
are any two points in the plane of the loop and B z is 
the operator for the magnetic field component normal to 
the loop. Both coordinates are integrated over the the 
surface of the loop to obtain the flux noise. In thermal 
equilibrium, the integrand can be written in the standard 
way in terms of the dyadic Green function and the mean 
thermal photon number nth- 

(B,(x,z;w)B^(y,z;w , )> = *(w-w , )^ 

7T 



xlm 



V x G(x,z;y,z;w) x V (h th + 1) . (2) 



If atoms are trapped near the superconductor with their 
spins parallel to the surface, this same noise in the mag- 
netic field component B z can drive spin flip transitions. 
For an atom at position ta, the rate is given by 
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(3) 

where is the magnetic dipole transition matrix el- 
ement between the initial and final Zeeman sublevels. 
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Here oj is the resonant transition frequency, correspond- 
ing to the atomic level splitting. The total spin-flip rate 
r is related to T z by T = i/2Y z . 

Neglecting the vacuum contribution to the spin flip 
rate and using the Weyl expansion for the scattering part 
of the Green function [HI, |3 , we find that 



V x G(x, z; y, z; w) x V 

d 2 h 
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where the perpendicular and parallel wavevector compo- 
nents are related by k z = \J U}2 1 c 2 ~ k 2 ancl r s is the 

Fresnel reflection coefficient for s-polarized (TE) waves, 
whose electric vector is perpendicular to the plane of in- 
cidence. 

In order to integrate over the circular pick-up loop of 
radius R, we use the identity 



. Thus, inserting Eq. J7]), together with 
the power series expansion ([S]), into Eq. ([T]), we obtain 
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The final approximation consists of assuming that the 
spin flip rate T z follows a strict power law with respect 
to the atom-surface distance: T z oc z~ n . This is cer- 
tainly the case in various limiting regimes [ill, El) H3| 
when the length scales relevant to the problem (tran- 
sition wavelength, atom-surface distance, skin depth of 
the substrate material etc.) can be well-separated. The 
derivatives in Eq. ((9]) then become 
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In this way, we perform the summation over s in Eq. Q 
to obtain 



where k|| ■ (x — y) = kulcosip with I = |x — y| and Jo(fciiZ) 
is the zeroth-order Bessel function. For the radial inte- 
gration we have 



I 



dllJ (k\ { l) = —J x {k\\R) , 
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giving the result 



d 2 xd 2 ye lk r(*-y) = A 2 



(k\\R) : 



:J?(hR) (7) 



for double integration over the pickup loop. Here we 
explicitly pull out the factor A 2 = (ttR 2 ) 2 , which is the 
squared area of the loop. 

The Weyl expansion ([4| of the Green function also re- 
quires us to integrate over transverse wave vectors. Since 
this cannot be done in closed form, we express the right 
hand side of Eq. ([7]) as a power series in ku Qjj], 
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In cases of practical interest, the distance z between the 
trapped atoms and the surface (typically l-100^m) is 
very small in comparison with the free-space wavelength 
of the spin-flip transition (typically 3cm-300m) . As a re- 
sult, the integral over ku is entirely dominated by the 
region in which k 2 3> u> 2 /c 2 , where k z « ik\\. When 
Eq. (@| is integrated to obtain the flux, in accordance 
with Eq. ([I]), the powers k 2s arising from the expansion 
© of the Bessel function can be obtained by differentiat- 
ing with respect to the atom-surface distance z, that is, 



S*(f) = 
X3F2 



h 2 A 2 
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where 3-F2 denotes a hypergeometric function. Equa- 
tion pip is the result we were seeking, connecting the 
measured flux noise spectrum S<p(f) to the anticipated 
atomic spin flip rate |T Z . The dependence of this con- 
nection on distance is controlled by the argument (R/z) 2 
of the hypergeometric function and by the power law as- 
sociated with the spin flip rate. 

In the limit of small R/z, when the size of the pick-up 
loop is small compared to its distance from the surface, 
the hy perg eometric function in Eq. pd|) can be approxi- 
mated [151 ] by 



3 r 2 
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Hence the flux noise spectrum is essentially proportional 
to the spin flip rate times the squared area of the pick-up 
loop. This reflects the fact that the loop is small com- 
pared with the transverse correlation length (of order z) 
and therefore the flux directly samples the local magnetic 
field. 

The flux measurements reported in [Llj were made in 
the opposite limit, R^> z, with a large pick-up loop lo- 
cated very close to the superconducting surface. Assum- 
ing a power law T z oc 1/z 4 (corresponding to the limit 
J « z with 6 bein g th e penetration depth of the sub- 
strate material [HI, Hi, [H, HJ), we obtain the limiting 
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FIG. 1: Reproduced from Ref. [14j|. Spectral density of mag- 
netic flux noise, S<s>(f), versus frequency / = uj/2tt for 15 
temperatures above the Kosterlitz-Thouless-Berezinskii tran- 
sition. Scatter at higher temperatures is due to subtraction of 
SQUID noise. Dashed lines have slope —1 and 0. Inset shows 
dV/dl versus T. 



R/z > 2.3, giving a spin flip rate of T z > x ■ 2 x 10 6 s _1 . 
At spin-flip frequencies above 10kHz (corresponding to 
a quantisation magnetic field stronger than 1/iT), the 
value of x given in Ref. fl4j is below 10~ 9 , corresponding 
to a trap lifetime in excess of 500 s. This rather slow T\ 
relaxation rate is very promising from the point of view 
of keeping atoms trapped near a superconducting sur- 
face. At the same time it may well be fast enough to be 
measured in the very benign environment of a cryostat. 

As well as inducing atomic spin flips, the magnetic 
field fluctuations can generate noise in the relative phase 
between Zeeman sublevels. In the presence of a static 
field Bo normal to the surface, the variance of the phase 
between levels 1 and 2 after time T is given by 



Mil) 2 



dt / dt' (B z (t)B z (t')) , (16) 



where B z {t) is the noise field and does not include the 
constant field Bq. This phase noise can be related to the 
spin flip rate. We connect B z (t) to B z {u) through 
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and we use Eqs. (J2]) and ([3]) to obtain 
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Comparing Eqs. p^)) and (fl"2")l we see that the spectral 
density of the flux noise measured in a large loop is also 
proportional to A 2 T Z but is suppressed by an additional 
factor of order (z/R) 3 . 

For the purpose of quantitative comparison, we take 
the magnetic dipole matrix element to be /U12 = Mb/2, 
corresponding to a transition between Zeeman sublevels 
\i) = \F = 2,m F = 2) and |/) = \F = 2,m F = 1) of a 
ground-state rubidium atom. Equation (jTTJ) then gives 



16 mi 



T z A 2 f(z,R), 



(14) 



where e and m e are the charge and mass of the electron 
and f(z,R) denotes the function in Eq. (|12| or Eq. (fT3|) . 
In order to make contact with the flux measurements 
reported in Ref. [l4| and reproduced in Fig. [TJ we write 
S<p(f) = x ■ $q, where $o = h/(2e) is the flux quantum. 
This gives 



rf 2 {B z (u,)Bt(u,')) = -r z (u;)6(u;~c,>) 



Substitution of Eqs. ([T7J) and {T 
yields the result 



into Eq. ([T 



(18) 
then 
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Now, T z (uj) is proportional to the measured flux noise 
S<j>(f) [Eq- (|15|]. which we know is constant up to a char- 
acteristic frequency /e, as shown in Fig.[T] Let us call this 
low frequency rate T(0). Above that, T z (ui) drops off as 
roughly 1/oj. The integral in Eq. (|19p is completely domi- 
nated by the low frequency range between < lu < 2tt/T, 
so we can make the approximation T z (u)) = T(0) pro- 
vided the observation time satisfies T > 1//$. Since 
exceeds 100 Hz, this is the case for T > 10 ms. Indeed, 
for large enough observation times the integral kernel ap- 
proximates the 8 function, 



2 (1 - coswT) 
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Equation (|19|) then reads simply 
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a[m]2 = (M22-M11) 2 Tz{Q)T 



(21) 



The SQUID loop that measured the flux in [14[ had an 
effective area of A = 180 x 900^m 2 w 2 • 10 _7 m 2 and had 



Supposing once again that states 1 and 2 are the 
ground states \F = 2,m F — 2) and \F — 2,m F = 1) of 
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a rubidium atom, the ratio of magnetic matrix elements 
squared has the value of unity and we obtain the partic- 
ularly simple result a[4>(T)] 2 = TjO)T ~ 2 x 10 6 x(0)T. 
Since the value of x(0) reported in [l4| is in the range 
10~ 9 — 10~ 5 , the corresponding dephasing lifetime T2 is in 
the range 50 ms — 500 s. This provides a very convenient 
time scale for the study of vortex noise using Ramsey in- 
terferometry, in which atoms prepared in a coherent su- 
perposition of Zeeman sublevels are later interrogated to 
measure the time-evolution of the coherence. The vortex 
field noise would be manifest as a loss of Ramsey fringe 
visibility with time, which could be studied as a function 
of the atom-surface distance. It should also be possible to 
explore the transverse coherence of the noise by varying 
the transverse extent of the cloud and measuring the loss 
of Ramsey fringe visibility as the cloud length increases. 

Measurements using cold atoms may also be able to 
image the vortices. The typical vortex separation of 
£ ~ 2/im 14], could be resolved by bringing the atom 
cloud to a similar distance from the surface, where the 
field of each vortex is of order <!>□/ (7r^ 2 ) ~ 1 G. At this 
close approach, the spin-flip lifetime is strongly reduced, 
but lifetimes approaching 100 ms may nevertheless be 
achieved. One imaging method would be to study the 
density distribution of the atoms, which is altered by the 
presence of the vortices through the effect of the vortex 



fields on the trapping potential. This approach is used to 
image classical current distributions in wires on an atom 
chip 0. The motion of the vortices could be tracked 
through the motion of the density patterns in the atom 
cloud. 

To conclude, we have shown that it is feasible to de- 
tect vortex dynamics in two-dimensional superconduct- 
ing films by means of trapped cold neutral atoms. In par- 
ticular we have considered the rate of atomic spin flips 
due to the magnetic field noise from the vortex motion. 
At 100 fJLTn from the surface we find that this lifetime can 
be in excess of 500 s, giving ample time to study the vor- 
tices. We have also considered the dephasing time for 
superpositions of Zeeman sublevels and find that this is 
short enough to be a sensitive measure of the vortex field 
noise. Finally, we have noted that spatial imaging of 
the vortices should be possible using cold atoms trapped 
close enough to the surface. These estimates show that 
cold atom clouds offer a sensitive new probe for the study 
of vortex dynamics in superconducting thin films. 
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